logo

Introduction

The U.S. Census Bureau collects and distributes data under a handful of different programs. Two of the more commonly used programs are the Decennial Census and the American Community Survey (ACS). The Decennial Census is a definite source of demographic data but only is collected every ten years; it includes a limited number of variables such as number of households and total population. The ACS is a program that provides data estimates on a one, three, and five year timeline; ACS data is collected more frequently but the data estimates have a margin of error that must considered because the data is taken from a small sample of the total population. The ACS includes many more variables compared to the Decennial Census that relate to transportation, income, and housing. Both the Decennial and ACS datasets have similiar data structures. Each row in both datasets include a particular variable and a number the indicates the total number of households or persons that characterize that variable.

Summary of TRPA Census Data

Population

  nv_acs <- get_acs(geography = "tract", year=2016, 
                variables =  c("B01003_001"),
                state = "NV", county=c("Washoe", "Douglas")) %>%
  mutate(data_source="2016 ACS 2016 5-year Estimate")
ca_acs <- get_acs(geography = "tract", year=2016, 
                variables =  c("B01003_001"),
                state = "CA",county=c("El Dorado", "Placer")) %>%
  mutate(data_source="2016 ACS 2016 5-year Estimate")
ca_decen <- get_decennial(geography="tract", variables= c("P001001"), 
                          state= "CA", year= 2010, county=c("El Dorado", "Placer")) %>% 
  rename(estimate=value) %>% mutate(moe=0, data_source="2010 Decennial Census")
nv_decen <- get_decennial(geography="tract", variables= c("P001001"), 
                          state= "NV", year= 2010, county=c("Washoe", "Douglas")) %>% 
  rename(estimate=value) %>% mutate(moe=0, data_source="2010 Decennial Census")
all<- bind_rows(nv_decen,ca_decen, ca_acs, nv_acs) %>%
  left_join(data.frame(tract), by="GEOID") %>%
    filter(!is.na(STATEFP)) %>%
  group_by(variable, data_source) %>% summarise(total=sum(estimate), moe=sum(moe)) %>%
  mutate(variable_name = case_when ( variable== "P001001" ~ "Total Population",
                                    variable== "B01003_001" ~ "Total Population"),
         total=format(total, big.mark=",", scientific=FALSE),
         moe=format(moe, big.mark=",", scientific=FALSE)) %>%
  select(variable_name, variable, total,moe, data_source)
datatable(all, extensions = 'Buttons',
rownames=F,options=list(dom='t', 
          columnDefs = list(list(className = 'dt-center', targets = 0:1))), 
  class = 'cell-border stripe', 
colnames = c('Variable Name', 'Code', 'Total', 'Margin of Error','Data Source'))

Housing

Income

income <- c(`Less than $10,000`= "B19001_002",
                    `$10,000 to $14,999`= "B19001_003",
                    `$15,000 to $19,999` = "B19001_004",
                    `$20,000 to $24,999` = "B19001_005",
                    `$25,000 to $29,999` = "B19001_006",
                    `$30,000 to $34,999` = "B19001_007",
                    `$35,000 to $39,999` = "B19001_008", 
                    `$40,000 to $44,999` = "B19001_009",
                    `$45,000 to $49,999` = "B19001_010",
                    `$50,000 to $59,999` = "B19001_011", 
                    `$60,000 to $74,999` = "B19001_012",
                    `$75,000 to $99,999` = "B19001_013",
                    `$100,000 to $124,999` = "B19001_014", 
                    `$125,000 to $149,999` = "B19001_015",
                    `$150,000 to $199,999` = "B19001_016", 
                    `$200,000 or more` = "B19001_017")
  nv_acs <- get_acs(geography = "tract", year=2016, 
                variables =  income,
                state = "NV", county=c("Washoe", "Douglas"), summary_var = "B19001_001") %>%
  mutate(data_source="2016 ACS 2016 5-year Estimate")
ca_acs <- get_acs(geography = "tract", year=2016, 
                variables =  income,
                state = "CA",county=c("El Dorado", "Placer"), summary_var = "B19001_001") %>%
  mutate(data_source="2016 ACS 2016 5-year Estimate")
income1<-rownames_to_column(data.frame(income), "name") %>% select(name) %>% pull()
all<- bind_rows( ca_acs, nv_acs) %>%
  left_join(data.frame(tract), by="GEOID") %>%
    filter(!is.na(STATEFP)) %>%
  group_by(variable, data_source) %>% 
  summarise(total=sum(estimate), moe=sum(moe),total_estimate=sum(summary_est)) %>%
  mutate(total=format(total, big.mark=",", scientific=FALSE),
         moe=format(moe, big.mark=",", scientific=FALSE),
         variable_name="Household Income") %>%
  arrange(factor(variable, levels=income1)) %>%
  select(variable_name, variable, total,moe, total_estimate, data_source)
datatable(all, extensions = 'Buttons',
rownames=F,options=list(pageLength = 16, dom = 'Bfrtip',buttons = c('csv','pdf'),
          columnDefs = list(list(className = 'dt-center', targets = 0:1))), 
  class = 'cell-border stripe', 
colnames = c('Variable Name', 'Definition', 'Number of Households', 'Margin of Error','Total Households' ,'Data Source'))

School Enrollment

Travel Mode to Work

Household Size

Demographics: Race

Browse all ACS Variables

Search through the list below to determine which variable(s) you want to analyze. You can download all of the variables

Browse all Decennial Variables

Search through the list below to determine which variable(s) you want to analyze.

TRPA Census Tracts

Download the Data Tahoe Open Data

TRPA Block Groups

Download the Data Tahoe Open Data